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1. INTRODUCTION 


MIDTRAN is a computer code which calculates transmittance and radiation for 
paths through the earth's atmosphere in the 2 - 5 um spectral region, The code is 
based on a marriage of the AFGL codes LOWTRAN3(1) and HITRAN(?) The overall 
structure and formats of LOWTRANS3 have been retained for the input parameters, 


path geometry and continuum transmittance components. The HITRAN data ‘ape is 
used to precalculate the spectral absorption properties of atmospheric molecules, 
which are stored on tape and then used as input data by MIDTRAN. The use of these 
data stored on external library tapes allows MIDTRAN to perform calculations for 
multilayered atmospheres at reasonable speeds and with a spectral resolution of 
0.05 cm! or better. Because of the library tape's present structure, MIDTRAN 
is best suited for paths which fall below 15 km altitude at some point, The code is 
written in FORTRAN and is compatible with CDC and IBM formats. For radiation 
calculations, the user has the option of specifying a background blackbody source 
of arbitrary temperature and emissivity and then calculating the radiation as seen 
through an atmosphere. In following the LOWTRANS structure, the user has the 
choice of six model atmospheres or radiosonde data and of different atmospheric 
paths, horizontal, vertical, or downward. Thus, MIDTRAN has flexibility for 
those systems studies which require atmospheric transmittance at lower altitudes 


while maintaining good spectral resolution. 


Section 2 gives a description of the MIDTRAN software and describes te 
MIDTRAN library tapes. Instructions for using MIDTRAN are in Section 3 and 


=F, Selby and R. McClatchy, ''Atmospheric Transmittance From 0.35 to 28.5 um; 
Computer Code LOWTRAN3,'' Report No, AFCRL-TR-75-0255, AFGL/OPI, 
Hanscom AFB, Mass., May 1975. 


2p, McClatchy, et al., ''AFCRL Atmospheric Absorption Line Parameters 
Compilation, '' AFCRL-TR-73-0096, AFGL/OPI, Hansom AFB, Mass., 
January 1973. 


comparisons to transmittance data are in Section 4. The appendices contain 
listings of MIDTRAN, a sample run, a list of variables, a MIDTRAN flowchart, and 
a listing of MRDAT, the program which generates the library tapes. 


The support furnished by the Naval Weapons Center and the Air Force Avionics 
Laboratory via the Air Force Geophysics Laboratories is greatfully acknowledged. 
The contract monitors are Mr. S. Ted Smith (NWC), Dr. R. Sanderson (AFAL), and 
Mr. B. Sandford (AFGL). Previous support of the Defense Advanced Research 
Projects Agency for the development of a preliminary version (MRDA) of MIDTRAN 
for use on an HP-2100 minicomputer?) 
sedes the earlier MRDA code. 


is acknowledged. The present code super- 


The MIDTRAN code can be directly used on either a CDC6600 computer or a 
minicomputer with virtual memory capability. Execution times on a CDC6600 for 
both transmittance and radiation is approximately 1.8 sec/wavenumber for a path 
transversing 11 model atmosphere layers and in steps of 0.01 ea (100 calculations / 
wavenumber). Calculational times on the PRIME 400 minicomputer are about 4 times 
slower. Total times for a calculation depend on the machine's tape read speed; 
considerable time can be spent by the PRIME in skipping over files to get to a 
spectral region near the end of a library tape. The code is still in the developmental 
stage. Qualified requestors may obtain copies of the code and library tapes from 


Aerodyne Research, Inc.; a charge will be made for tape duplication. 


3p, Kryger and D. Robertson, ''MRDA - A Medium Resolution Data Analysis Code 
for the HP2100 Minicomputer, 't AFGL-TR-77-0044, AFGL/OPR Hanscom AFB, 
Massachusetts 01731. 


2, DESCRIPTION OF MIDTRAN SOFTWARE 


2.1 The MIDTRAN Code 


MIDTRAN is designed to make medium resolution atmospheric transmission 
and radiation calculations in the range of 1800 to 6000 cem7! over a wide variety of 
geometrical paths. The first part of the calculation consists of predicting the con- 
tinuum transmissions (H,O, No, molecular scattering) along with aerosol _absorp- 
tion for the chosen path and wavenumber interval. These calculations are carried 
out by what is essentially LOWTRAN3") with all the spectral calculations removed. 
In the second part of MIDTRAN, the medium resolution spectral calculations are 
performed, Magnetic tapes which contain a complete library of extinction coefficients 
are used to provide the data for computing the spectral contributions. After having 
calculated the totaltransmission due tothe spectral structure, the program then 
combines the continuum and medium resolution results; if desired, the radiation 
is also calculated at each frequency. The frequencies, radiances, and transmittances 
are then written to a disk file, associated with FORTRAN logic unit 9. In the third 
part of MIDTRAN, this disk file is rewound, and the transmittances and radiances 
are then degraded to the desired spectral resolution using the available slit function 
before being printed out and/or plotted. 


In the process of computing the continuum results, intermediate values are 
saved for later use in the spectral calculations, For example, the pressure, 
temperature, and altitude of each layer traversed by the geometric path are stored. 
In addition, the transmission through each layer is also stored in an array TRANI for 
use in calculating the radiation. Finally, the atmospheric concentration of HO and 
Og in each layer along with the molecular density (of all gases) for the particular 
path through each layer is saved. Using this information from the continuum part of 
the calculation and data from the library tape, MIDTRAN calculates the spectral 


transmission and radiation over the geometric path at frequencies defined by the 


Spauiwneestaeee 


input, The spectral absorption coefficients are read from the library tape. They 


are tabulated for the six important infrared atmospheric molecules: 


H,O, CO O., N,0, CO, and CH,. 


9? 
The tapes are organized in 2 cm! blocks over the 1800 to 6000 em~! range 

that MIDTRAN operates. In each wavenumber block, the extinction coefficients for the 

six molecular species are tabulated at 9 pressure-temperature points, the wave- 

numbers being chosen to represent the structure of the absorption spectra for the 

particular species in that wavenumber block. For each species, the wavenumbers 

were selected so as to define the spectra by identifying the most important lines in 

the region. In combining the separate contributions from far and near spectral 

lines, the extinction coefficients were calculated near the line center at 10 points, 

spaced 0.01 em7! apart, andat 0.1 cm-1 irtervals between adjacent strong lines. 

In order to obtain the extinction coefficient at a particular pressure, temperature 

and frequency, the program performs linear interpolations over the pressure/ 


temperature matrix and then over frequency. 


Presently, the slit function library in MIDTRAN contains the option for no 
slit function at all and for a generalized slit function which requires the user to 
input two arrays for its definition, the slit width and the shift. The plotting option 
requires that the user input certain titles, initial axis values, and scaling 
parameters for use by the plotting software. The plotting software now in the 
program is designed for a PRIME 400 system with a Versatec printer/plotter. 
The user must examine this part of the code to determine its compatability with 
his system. Printed output is columnar and is blocked in sections which contain 


a maximum of 240 pairs of output data. 


2.2 The MIDTRAN Library Tapes 


The MIDTRAN Library tapes contain the spectral absorption coefficients for the 


six atmospheric species which have significant absorption in the 1800 - 6000 em7l 


region, The species are: H,O, CO4,0., N,0, CO, and CH,. The absorption 


Waa ee 


coefficients are calculated at selected (P, T, v) points and then written onto a tape 
that is accessed by MIDTRAN. The CDC6600 computer at AFGL was used to 


generate these tapes. 


2,2.1 Choice of Spectral Absorption Coefficients 


The data in the MRDA library tape are organized so as to define the absorption 
spectra for the species in as compact a form as possible, Thus, nine pressure- 
temperature (P, T) points are used to describe the atmosphere, and the total 


number of wavenumber points within each block is limited to 250, 


The choice of the (P, T) points is based on the expected range of atmospheric 
pressures and temperatures. Figure 1 shows the (P, T) variability, along with 
(4) The 
heavy dots within the circles show the nine (P, T) points at which the spectral 


illustrative radiosonde data, taken from several AFGL Mission. 


absorption coefficients are calculated. Pressure/temperature points for pressures 
below 100 mb are not included in the tape at this time, since the dominant part of 
the atmosphere is at lower altitudes. The tape program contains Doppler line- 


shapes, so the user can generate his own high altitude tape to use with MIDTRAN. 


A considerable savings in the total number of wavenumber points at which the 
spectral absorption coefficients must be stored is obtained by identifving the 
stronger spectral lines within each wavenumber block. When one or more of the 
species have strong lines within a block, the absorption coefficients are calculated 
at the peak, at 10 points about line center 0.01 em] apart, and at intervals of 
0.1 em7! between adjacent peaks. Spectral absorption coefficients for intermediate 


values are obtained by linear interpolation, 


“a. Sandford, et al., ''Aircraft Signatures in the Infrared 1,2 to 5.5 Micron 


Region, '' AFGL-TR-76-0133, Air Force Geophysics Laboratory (OPR), Hanscom 
AFB, Mass., 01731, (June 1976). 
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Figure 1 - Temperature and Pressure Variations of the Atmosphere. 
Radiosonde Data from Several AFGL Measurements are Indicated. 
The Outer Lines Indicate the Approximate Range of Atmospheric 
Temperature Fluctuations, The Center Line is the U.S. Standard 
Atmosphere, 
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2.2.2 The MRDAT Code 


The MRDAT code is operational on the AFGL CDC6600 computer. It calculates 
the spectral absorption coefficients from line parameters contained on the HITRAN 


(2) and writes them to an external tape. These tapes comprise the MIDTRAN 


data tape 
tape library. The Lorentz lineshape function is used to describe the contribution of 
overlapping line tails. In the 2360 - 2500 em7} spectral region, the co, lineshape 
includes the Burch form factor) as modeled by Kaplan, et al (6) 4 Voigt 

lineshape is included in the program; internal logic selects this lineshape when the 


Doppler width becomes comparable to the Lorentz width. 


The absorption coefficients are calculated in two steps. The contribution of 
lines external to each wavenumber block are calculated at two points (the edges of 
the block); linear interpolation is used to determine their contribution at intermediate 
frequencies, The contribution of the nearby lines is calculated at each wave- 
number point within the block. The two results are combined and then written onto 


the library tape. 


The spectral absorption coefficients for each molecule are calculated for 
1 cc of pure gas at STP (i.e., 2.69 x es molecules). MIDTRAN includes the 
concentration when calculating the transmittance. The codes are written for 
9 (P, T) points and 6 atmospheric molecules, They are specified in the MRDAT 
input cards. So long as these parameters maintain the structure shown in Fig. 1, new 
tapes which are tailored to specific problems (like high altitude) are easily generated. 


2.2.3 MRDAT Input Parameters 


The spectral absorption coefficients are calculated from the molecular line 
parameters on the AFGL HITRAN data tape. A listing of the program (MRDAT) is 


"DE. Burch, D.A. Gryvnak, R.R. Patty, andC.E, Bartky, Journal, Optical 
Society of America, 59, 267, 1969, 


2D, Kaplan, M.T. Chahine, J. Susskind, J.E. Searl, Applied Optics, 16, 
322, 1977. 
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given in Appendix B. The data for the input parameters are given by the following 


sequence of read input lists: 


Ks NPTPTS, MSPEC Format (812) 
2. P(), l=1, NPTPTS Format (8E10,.0) 
3. T (I, I1=1, NPTPTS Format (8E10.0) 


4. W(M, M- 
5. V1, V2,DV, VLWST, VHGHST, DELTV, BOUND Format (7E10.3) 


a ey Format (7E10.3) 


6. SSTR, VBLOCK, DV2 Format (3E10, 3) 


The input quantities are: 


NPTPTS 
32: 
T 
Ww 


Vi 

V2 

DV 
VLWST 
VHGHST 
DELTV 


BOUND 


SSTR 
MSPEC 


number of (P, T) points 

pressure values 

temperature values 

species column density = 0.269E20 molecules/em” for 
(H,O, CO,, O., N,0, ECO, CH, O,)* 
lower frequency limit (em71) of the library tape 

upper frequency limit (em7}) of the library tape 
frequency increment for calculating between strong lines 


lower frequency bound, em 1, for consideration of distant 
lines, (presently, overidden internally) 


upper frequency bound, em=1, for consideration of distant 
lines, (presently, overidden internally) 


frequency increment for distinguishing between near and far 
lines, 


distance (em=1) from line center beyond which a line is not 
included, presently fixed at 20.0 cm 


lower line intensity limit used for accepting lines. 
identification of the six molecules (Set = 123456) (2) 


* 

Since oxygen does not have any important absorption bands below 6000 em=i, it is 
not included as one of the six species in MIDTRAN, but could be used in place of a 
molecule on a new tape library. 


VBLOCK = frequency interval length (em"}) into which the range [vi, v2] 
is divided for blocking. 
DV2 = frequency increment (usually 0.1 em7!) for calculations 


between strong lines. 


Figure 2 is a schematic which illustrates the choice of these wavenumber 


i. 


parameters. 


AL-78-303 


| VE LOCK 


DELTV DELTV 


''Far'' Lines ate ‘Near’! Lines ————- Wart! Lines 


Figure 2, Schematic Showing the Definition of Wavenumber Parameters 
for MRDAT,. Note thet the Second Line is Weaker Than SSTR 
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3. OPERATING INSTRUCTIONS 


3.1 Input Data and Formats 


Many of the input cards and operations are very similar to those of LOWTRAN. 
In particular, the card input sequence is nearly identical. The input cards can 
be divided into two blocks, specification of the case to be calculated and specification 
of slit function and/or plotting parameters and formats. Up to 10 cases can be 
run. All calculations (including multiple cases) are completed in the first block 
(CARDS 1 - 4) before reading the second block of data (CARDS 5 - 7). Cards 
1, 3, and 4 are identical to LOWTRAN3 input cards") 


The data necessary to specify a given problem are given in a four card 


sequence as follows: 


CARD 1 MODEL, IHAZE, ITYPE, LEN, JP, IM, M1, M2, M3, 
NLDAT, R_ Format (1013, F10.3) 
CARD 2 IRAD, EMIS, TBACK, NTS,NTP, NRS, NRP, XORG, YORG 


Format (10, 2F10.3, 415, 2F10.3) 
CARD(s) 2A (Atmosphere data cards when M = 0 or M = 7) 


CARD 3 H1, H2, ANGLE, RANGE, BETA, VIS 
Format (6F10.3) 
CARD 4 V1, V2, DV 
Format (3F10.3) 
CARD 1 (Model = -1 to indicate last calculation) 
CARD 5 TITLE 


Format (20A4) 


(Slit Function Parameters) 


CARD 6A WIDTH, SHIFT, NS 
Format (2F10.5, 110) 


CARD 6B XSS (I), I= 1, NS 
Format (8F10.5) 


CARD 6C SS (I) l= 1, NS 
Format (8F10.5) 


14 


(Plotting Parameters) 


CARD 7A XTITLE 
Format (20A4) 


CARD 7B YTITLE 
Format (20A4) 
CARD 7C XAXIS, XINIT, XSCALE, DXT, NMINX 


Format (4E10.4, I10) 


CARD 7D YAXIS, YINIT, YSCALE, DYT, NMINY 
Format (4E10.4, 110) 


If MODEL = 0 or 7, meteorological data used to describe the atmosphere are 
inputted on CARD(s) 2A. Transmittance and radiation calculations for all the various 
cases are completed with results written to an external file, before either the slit 
function is used or the plotting routine is employed. The external file is associated 
with FORTRAN logical unit 9. CARDS 1 - 4 can be repeated to perform up to ten 
calculations, ending with a MODEL = -1 on CARD 1. Another cyclical sequence 
of-input data follows this card to specify the title, slit function parameters, and 
plotting parameters for each of the cases. Up to two plots (radiation and transmittance) 
can be made for each case and card set. The first block (CARDS 1 - 4) is described 
in Subsection 3.2 and the second block in Subsection 3,3. 


3.2 Input Parameters 
3.2.1 CARD1: MODEL, IHAZE, ITYPE, LEN, JP, IM, M1, M2, M3, NLDAT, R 


The parameter, MODEL, selects one of the six geographical model atmospheres," 


specifies that meteorological data are to be used in place of the standard models, or 


indicates the end of the first blocks of data (the second block being the output parameters). 


THAZE specifies whether aerosol attenuation is to be included in the calculation 
or not. For any problem, the atmospheric path must be specified as one of three 
types according to ITYPE and LEN. The rest of the quantities given on CARD 1 (which 


can be left blank if not required) provide the user with options to suppress printing 


(JP), to intermix the six standard model atmospheres (M1, M2, M3), to input 
anew model atmosphere (IM, NLDAT), and to specify the earth radius (R). The 
options for the above parameters and their use are described below: 
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MODEL = -1 indicates end of first data blocks 


= 0 indicates meteorological data are specified for a horizontal 
(constant pressure) path. 


selects TROPICAL MODEL ATMOSPHERE | 
selects MIDLATITUDE SUMMER 
selects MIDLATITUDE WINTER ‘ 
selects SUBARCTIC SUMMER / 
selects SUBARCTIC WINTER | 
selects 1962 US STANDARD | 


indicates a new model atmosphere (or radiosonde data) is to i 
be inserted | 


iH} 
NAannrt WwW dw 


IHAZE = Q means no aerosol attenuation included in the calculations, 


= 1 or 2 if aerosol attenuation is required (see also, CARD 2). 


If IHAZE is set equal to 1 or 2 and visual range (VIS) is not specified on | 
CARD 2, the program will then automatically select visual ranges of 23 or 5 km, 
respectively. 


ITY PE = 1 for a horizontal (constant pressure) path. 
= 2 for a vertical or slant path between two altitudes. 


= 3 for a vertical or slant path to space, 


The TYPE 1 path should not be confused with a 90° path where the local height 
at the endofthetrajectory is significantly different from thht at the beginning, In 
such a case, specify the path according to ITYPE = 2. 


LEN 


0 for normal operation of program. 


= 1 selects the downward TYPE 2 path shown in Figure 3(e). 


The parameter LEN, can be ignored (that is, left blank) for the majority of 
cases. It need only be used for a downward looking path (H2< H1) when two paths 


are possible for the same input parameters. In such a case, a computer printout . 
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fa 


Hi 


(a) (b) (c) 


H MIN 


SHORT PATH LONG PATH 


(d) (e) 


Figure 3, Geometrical Path Configuration for: (a) Horizontal Paths (Type 1), 
(b) Slant Paths Between Two Altitudes H1 and H2 (Type 2), and 
(c) Slant Paths to Space (Type 3). For Downward Looking Paths 
Where HMIN < H2 < Hl, Two Trajectors Are Possible As 
Indicated in (d) and (e). The Angle @ Corresponds to ANGLE on 
CARD 3. From Ref, (1). 
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eee 


statement will be given, indicating that the user has two choices for the problem 


and the shorter path (see Figure 3(d)) has been executed. Set LEN = 1 for the 


longer path case. 


JP 


IM 


NLDAT 


il} 


on FF © 


for normal operation of program. 


1 printout of data. 


additional printout, including a 0.1 cm™ 
partial printout at a resolution of DV. 


highest level printout includes absorption coefficients from 
tape. 
when radiosonde data are to be read in initially. 


for normal operation of program or when subsequent 
calculations are to be run with MODEL = 7, 


number of levels to be read in for MODEL = 7, 


Note that IM and NLDAT are only used when MODEL = 7 and then only on 
the first calculation when the atmospheric data are read from Card(s) 2A. 


The parameters M1, M2, and M3 can each take any integral value between 0 and 
6. Se: M1 = M2 
supplement the altitude profiles of temperature, water vapor, and ozone, respectively, 


for any given atmospheric model specified by MODEL. For example: 


M1 


M3 


M3 = 0 for normal operation of program, They modify or 


selects the TROPICAL temperature altitude profile 


selects the MIDLATITUDE SUMMER temperature altitude 
profile 


selects the 1962 US STANDARD temperature altitude profile 
selects the TROPICAL water vapor altitude profile 


selects the MIDLATITUDE SUMMER water vapor altitude 
profile 


selects the 1962 US STANDARD water vapor altitude profile 
selects the TROPICAL ozone aititude profile 
selects the MIDLATITUDE SUMMER ozone altitude profile 


Sear erat ENTER nara tegen en nas 
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: 


padi ms = 


R = radius of the earch (km) at the particular geographical 
location at which the calculation is to be performed. 
If R is left blank, the program will use the midlatitude value of 6371.23 km 
when MODEL is set equal to 0 to 7. Otherwise, the earth radius for the appropriate 
standard model atmosphere (specified by MODEL) will be used. 


When MODEL = 0 or 7, the new atmosphere (model or radiosonde data) 
is inserted between CARDS 2 and 3. 


This card is in addition to the LOWTRAN input cards. It determines whether 


3.2.2 CARD 2: IRAD, EMIS, TBACK, NTS, NTP, NRS, NRP, XOR, YOR 
atmospheric radiation is calculated and specifies the slit function and plot 


parameters. 
IRAD = 1/0 if radiation calculations are/are not to be made. 
EMIS = emissivity of a background radiation source located at 

| the beginning of the path. (Include, if IRAD= 1.) 

TBACK = temperature (in degrees Kelvin) of the background 
radiation source, (Include, if IRAD= 1.) 

NTS = 0 when using a previously specified value. 

= +/-1 use variable slit function on transmittance and plot/don't 
plot. 

, = +/-3 don't use any slit function (leave points as they are) for 
transmittance and plot/don't plot. 

NTP = plot transmittance vs cm-l, 

| 

=. plot transmittance vs microns. 

| NRS = 0 +1, +3 (same as NTS, for radiation). 

NRP = +] (same as NTP, for radiation). 

| XOR = coordinates (in inches) for the lower left corner of the 

| YOR plot. 
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The background radiation source is calculated using the temperature dependence 
of the blackbody function and a surface emissivity given by EMIS. Radiation from this 
gray body source is then propagated through the atmosphere from H1 to H2. 


It should be noted that all transmittance and radiation calculations for all of 
the various cases are made and written to an output file, before any slit functions 
are used or plotting software is employed. The external file is associated with 
FORTRAN logical unit 9, After the transmittance and radiation calculations for all_ 
of the cases have been written to unit 9, it is rewound and used as input for the slit 
function and plotting subroutines, For degrading and plotting, the order of processing 
then proceeds as follows: CASE 1 transmittance, CASE 2 transmittance, ...., 
CASE N transmittance, CASE 1 radiation, CASE 2 radiation, ...... , Case N 


radiation, This is illustrated in the following matrix: 


NTS (CASE 1), NTP (CASE 1) 
NTS (CASE 2), NTP (CASE 2) 

z Briefly, Column 1 monitors 
NTS (CASE N), NTP (CASE N) the present slit function being 
NRS (CASE 1), NRP (CASE 1) used and whether or not to 

plot. Column 2 determines 
NRS (CASE 2), NRP (CASE 2) the units if the user has 


chosen to plot. 


NRS (CASE N), NRP (CASE N) 


Once a particular slit function has been specified in Column 1; under the NTS 
or NRS parameter, the NTS/NRS column can be left blank until a new slit function 
is to be used (with the exception that in going from plotting to no-plotting or vice versa 


the NTS/NRS parameter has be to explicitly entered). 


20 


Likewise, the NTP/NRP parameter can be left blank after being specified, 


until a new set of plotting units is desired. However, the NTP/NRP parameter need 


not be respecified following a series of no-plot options under the NTS/NRS parameter. 


3.2.3 CARD2A: (For MODEL = 0 or 7) 


If MODEL = 0 and ITYPE = 1, then meteorological data for a horizontal 
(constant pressure) path are to be inserted between CARD 2 and CARD 3 
as follows: 
H1,P, T, DP, RH, WH, WO, VIS, RANGE 
Format (3F10.3, 2F5.1, 2E10.3, 2F10.3), 


where the above parameters refer to altitude (km), pressure (mb), ambient 
temperature (°C), dew point temperature @C), relative humidity (7), water vapor 
density (gm m73), ozone density (gm m7), visual range (km), and path length (km), 
respectively. It is only necessary to specify the quantities underlined with the solid 
line and one of the quantities underlined with the dashed line. The ozone density WO 
can be specified using the parameter M3 on CARD 1, if data are not available. In 
the latter case, a value will be calculated at altitude H1, based on the appropriate 


model atmosphere selected by M3, 


If MODEL = 7 and IM = 1, then a new model atmosphere must be inserted at 
this point, between CARD 2 and CARD 3 (1) The number of atmospheric levels 
to be inserted is given by NLDAT on CARD 1. The format for atmospheric data at 
each of the levels is: 


Z, P, T, DP, RH, WH, WO, AHAZE 


ad — 


Format (3F10.3, 2F5.1, 2E10.3, F10.3) 


The first level should be at Z= 0.0. These parameters are the same as defined 
above in this subsection, excepting AHAZE, the aerosol number density (em7}), 
It is only necessary to specify these quantities underlined with a full line and one 
of the quantities underlined with the dashed line. If the aerosol number density was 


not measured as a function of altitude and the user wishes to include aerosol 
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attenuation in the calculation, set IHAZE = lon CARD1. In this case, MIDTRAN 
will use the aerosol models already contained in the program and interpolate to give 
aerosol number density values at the same altitudes as the radiosonde (or new model 
atmosphere) data. The program will then look for a sea level visual range (VIS) to 
be specified on CARD 3. If VIS is not specified, a 23 km sea level visual range 

will be assumed. If aerosol attenuation is not required, set IHAZE = 0 on Card1 


as before. 


3.2.4 CARD3:H1, H2, ANGLE, RANGE, BETA, VIS 


CARD 3 is used to define the geometrical path parameters for a given problem. 


H1 = initial altitude (km) 

H2 = final altitude (km) 

ANGLE = initial zenith angle (degrees) as measured from H1 
RANGE = path length (km) 

BETA = earth center angle subtended by H1 and H2 (degrees) 
VIS = sea level visual range (km) 


It is not necessary to specify every quantity given above, only those that 
adequately describe the problem according to the parameter ITYPE (as described 


below). 


(1) Horizontal Paths (ITYPE = 1) 
(a) specify H1, RANGE, and VIS only 


(b) if nonstandard meteorological data are to be used (that is, is 
MODEL = 0 on CARD 1), then the radiosonde data must be 
specified on CARD 2A and CARD 3 is omitted, 


(2) Slant Paths to Space (ITYPE = 3) 
(a) specify H1, ANGLE, and VIS 


(b) specify H1, HMIN, and VIS (for limb viewing problem where HMIN 
is the tangent height or minimum altitude of the path.) 


(3) Slant Paths Between Two Altitudes (ITYPE = 2) 
(a) specify H1, H2, ANGLE, and VIS 
(b) specify Hl, ANGLE, RANGE, and VIS 
(c) specify H1, H2, RANGE, and VIS 
(d) specify H1, H2, BETA, and VIS 


For cases (b) and (c), the program will calculate H2 and ANGLE assuming 
no refraction and then proceed as for case (a). This method of defining the 
problem should be used when refraction effects are not important; for example, consider 
ranges of a few tens of km at zenith angles less than 80°. It can also be used for 
larger angles (including 90°) provided that the path lies within one atmospheric 
layer. 


Leave ANGLE and RANGE blank in case 3(d). This method can be used when the 
geometrical configuration of the source and receiver is known accurately, but the 
initial zenith angle is not known precisely due to atmospheric refraction effects. 


BETA is most frequently determined by the user from ground range information. 


In the cases of 2(b) and 3(d) above, the subroutine, ANGLE, is called in 
the program to determine the appropriate input zenith angle by the LOWTRAN3 
iterative technique ™) that takes atmospheric refraction into account. 


3.2.5 CARD4: V1, V2, DV Format (3F10.3) 


The spectral range over which transmittance data are required and the spectral 


increments at which the results are calculated is determined by this card. 


V1 = initial frequency in wavenumbers (em?) 
V2 = final frequency in wavenumbers (cm7}) where V2 > V1 
DV = frequency increment (or step size) (em74) 


Note that v = 10°/r where vy is the frequency in em7! and 2 is the 


wavelength in microns, 


For lower altitude paths, values of DV around 0.05 cm7! give sufficient 


accuracy; for high altitudes, 0.02 or even 0.01 em” is necessary. 


This completes the set of cards necessary to specify one transmittance/ 
radiation calculation, If more cases are desired, repeat the sequence. If no more 
cases are desired, CARD 1 with MODEL = -1 is inserted after CARD 4 and before 


the slit function/plot cards that are described in the next section. 


3.3 Output Parameters 


In the same way that the input data blocks are given in sequence, another 
sequence of data specifying the output format and parameters must be given, The 
MODEL = -1 card separates the two groups of cards. In the second sequence, the 
title, slit function, and plotting cards for each transmission or radiation calculation 
must be specified. Generally, one complete cycle in the sequence from this 


second section of input data is structured as follows: 


CARD 5 TITLE 
CARDS 6A - 6C SLIT FUNCTION PARAMETERS 
| CARDS 7A - 7D PLOTTING PARAMETERS 


\ 


All of these calculations are executed in one call to subroutine LIB, which is made 


just prior to stopping above statement #27. 


3.3.1 CARD 5: TITLE (20A4) 


The title is used on the plot and is printed with no change. 


3.3.2 CARDS 6: 


WIDTH, SHIFT, NS Format (2F10.5, 110) 
XSS(1), I= 1, NS Format (8F10.5) 
SS(I), I= 1, NS Format (8F10.5) 


WIDTH = width of slit function in em7} 


SHIFT = distance (em7?) between points at which the slit function is 
calculated 

NS = number of (XSS, SS) points to define the slit function (max. 8) 

XSS = wavenumber coordinate of slit function points 

ss = _ weighting function values for slit function 


When no slit function is desired (i.e., print the results directly), Cards 6 
are omitted, This is determined by setting NTS or NRS equal to +3 in Card 2, 
An illustrative example for the generalized slit function is given in Figure 4. The 
results are degraded to the desired resolution by integrating over the slit function. 
If the same slit function is used for subsequent calculations, Cards 6 are omitted, 
This omission must be reflected by zeros for the NTS and NRS parameters in Card 2. 


Arbitrary values of SS and XSS can be used, because the slit function is normalized. 


SS AL- 78-304 


XSS 
WIDTH 


Wavenumber (cm7!) 


Figure 4, Example of a Generalized Slit Function Specified By 
Six Values of XSS and SS 
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3.3.3 CARDS 7: 


XTITLE Format (20A4) 
: YTITLE Format (20A4) 
4 XAXIS, XINIT, XSCALE, DXT, NMINX Format (4E10.4, 110) 
YAXIS, YINIT, YSCALE, DYT, NMINY Format (4E10.4, 110) 
XTITLE = Title for abscissa units 
YTITLE = _ Title for ordinate units 
XAXIS = Length of x-axis in inches 
XINIT = Value of x at the origin 
SCALE = Change in value of x per inch of plot 
DXT =  x-units between major tic marks 
NMINX = Number of minor tic marks between the major ones. 


(Same definitions for y-axis. ) 


These plot parameters are those that are necessary to specify a plot on the 
PRIME 400 computer at Aerodyne Research, Inc. The user must modify these 
parameters and their definition in order to be compatible with his plotting software. 


The external plotting subroutines which MIDTRAN expects to find are: 


AXIS, 
PLOT, 
INIT, and 
ENDPLIT. 


They are called from subroutines LIB, FRAME, and PROUT. 
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4, COMPARISON TO DATA 


Comparison of atmospheric transmittance calculated using MIDTRAN is made 


to data taken by wR 


data's spectral structure. Figure 5 shows NRL data taken for a 5.12 km sea level 


The comparisons show that the code is able to calculate the 


horizontal path at a spectral resolution of 0.08 cm7! and a MIDTRAN calculation 
performed at a resolution of 0.01 em7! and degraded to 0.08 cm7! using a triangular 
slit function. The calculations combine the molecular spectral structure of a 
HITRAN calculation and the H9O and No continuum components from LOWTRAN. 


This yields a calculated spectra that compares favorably with the data. 


Figure 6 shows a comparison to the same data but on an expanded scale in 
the 2385 - 2450 cm7! spectral region. This illustrates the fall-off in the transmittance 
as one moves in towards the 4.3 ym CO9 band. The COs spectral absorption 
coefficients in this region include the tail contribution from all lines in the CO» band. 
The form factor of Kaplan, et ai ®) was used to modify the Lorentz lineshape. 
The calculations have the correct roll-off as exhibited by the data but underestimate 
the strength of this effect. Since the difference between the calculated and measured 
spectra decreases as one moves away from the CO, band, the difference is most 
likely due to the modeling of the CO» form factor. The parameterization used by 
Kaplan, et al. is already larger than the form factor calculated by Burch, et a) 
from their data, so no further changes were made to fit these data. Further 
theoretical studies and measurements at different atmospheric temperatures are 
required to better parameterize the lineshape in this region (and to verify that the 
strength is really due to COg and not to unexpected spectral structure in the No 


and/or H2O continuum components). 


TK, Haught and J, Dowling, ''Long Path High Resolution Field Measurements of 
Absolute Transmission in the 3.5 to 4.0 4m Atmospheric Window, "' Optics Letters, 
i, Eat (1977), 
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Figure 6, Comparison of NRL Data and MIDTRAN Calculation for a 
5.12 km Path at Sea Level From 2385 to 2450 em-1 


Figure 7 shows a comparison of AFGL data and MIDTRAN calculations for the 
transmittance to space from an altitude of 8.5 km. The transmittance data, which 
were taken from the AFGL KC135A flying laboratory, are obtained by measuring 


the solar spectrum and then dividing out the solar irradiance to obtain the atmospheric 


Local radiosonde data were obtained. The spectral resolution of the AFGL inter- 


(4) 


The MIDTRAN calculation was done at a spectral resolution of 0.01 em! 


then degraded to 3.8 em using the actual slit function of the AFGL interfereometer. 


The data were taken in the vicinity of Johnston Island in the Pacific. 
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in the calculation. The radiosonde is not able to measure atmospheric water vapor 


| 
1 
model atmosphere’ ) for the high altitudes. Table I shows the layers which were used | 


concentrations when the dew point depression temperature becomes less than -30°, 

Thus, the radiosonde inputs resulted in too much H,O absorption around 2000 em7l, 
the much lower H20O concentrations listed in Table I bring the calculations and data 
into better agreement. The ozone profile is that of the Tropical Atmospheric model. | 


TABLE I - Model Atmosphere Used for the MIDTRAN Calculation Shown in Figure 7. 
It is Based on Local Radiosonde Data Supplemented by the Tropical 


Model Atmosphere. 
Pressure Temperature 03 3 
(mb) (C) gm/m 
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APPENDIX B 
MRDAT PROGRAM LISTING 


PROGRAM MRDAT(INPUT,OUTPUT, TAPE2 , TAPE4 , TAPE6=OUTPUT, TAPES=INPUT, 


PROGRAM MRDAT (INPUT ,OUTPUT,TAPE2 , TAPE4 , TAPE6=OUTPUT, TAPE5=INPUT, 
X TAPE11,TAPF12,TAPE13) 


MAY 9 77 HITRAN MODIFIED FOR BLUE CO2 ARC TARLES 
JUNE 23 -- VOIGHT PROFILE AND FORM FACTOR 
JULY 77 MODIFIED FOR RED CO2 REGION 
AUG 77 REWRITTEN FO MRDA 
LOGICAL LCHK,LOGIC 
DIMENSION SUM1(6,9,1@1),SUM2(6,9,101) ,OMEGB(191), 
KSPECIE(7) .AMASS(7) ,JCALC(7) 
DIMENSION CAY3(6,9),CS1(9),CS2(7,9),CA(9) 
COMMON/INPUT/ P(18),T(18) ,W(7),V1,V2,DV,VLWST,VHGHST,DELTV, 
KBOUND,NPTPTS ,MSPEC(7) ,SSTR,VBLOCK ,DV2 
COMMON/OMG/ OMEGA(281,6),STOR(6,9,3084) 
COMMON/LINES/ GNU(25@),S(250) ,ALPHA(259) ,EDP(258) ,MOL(25@), 
KLCHK (258) ,TI(258,9) ,ITI(259) , TMAX 
COMMON/BLOCK1/MAX,CS1,CS2,CA,SUM1,SUM2 ,OMEGB,JMAX,NFILES 
DATA AMASS/18.,44.,48.,44.,28.,16.,32./ 
DATA SPFCIE/4HH20 ,4HCO2 ,4HO3 ,4HN20 ,4HCO ,4HCH4 ,4HO2 / 


PI=3.14159 
PP=1013. 
TO=2596. 


MAX=5 
VSTEP=59.8 
SLOWER=1.E8-27 


ROOT OF MOLECULAR WT FOR DOPPLER LINESHAPE 
DO 1 M=1,7 
AMASS (M) =SQRT(AMASS (M) ) 


THESE DEFINTIONS ARE NOT ACCURATE 
READ INPUT PARAMETERS (P=PRESSURE), (T=TEMPERATURE), 
W(1)=H28, W(2)=C@2,W(3)=83,W(4)=N28,W(5)=CO,W(6) =CH4 ,W(7)=82. 
Vi AND V2 ARE FREQUENCY LIMITS FOR WHICH OUTPUT RESULTS ARE REQUIR 
DV IS MONOCHROMATIC FREQUENCY INCREMENT. 
BOUND IS THE FREQUENCY FROM ANY LINE CENTER BEYOND WHICH THE LINE 


DEFINE MOLECULAR SPECIES 
READ INITIAL PARAMETERS 
CALL DATIN 
WRITE(4,3) V1,V2,NPTPTS 
FORMAT(2F1@.2,15) 
WRITE(4,4) (P(I),I=1,NPTPTS) 
WRITE(4,4)(T(I),I=1,NPTPTS) 
4 FORMAT(9F1@.2) 
*** CALCULATE CONTRIBUTION OF DISTANT LINES 
“axe PIRST DETERMINE OMEGBS (CALCULATIONAL POINTS) 
5 VBOT=FLOAT(INT(V1)) 
VTOP=FLOAT(INT(V2) ) 
IF(VTOP.GT.2360.8.AND.VBOT.LT.2360.8) VTOP=2360.8 
IF (VTOP.LT.V2) VTOP=VTOP+1.@ 
JMAX = INT( (VTOP-VBOT) /DELTV+1 .998@1) 
JMAX MUST BE .LE. 101 
IF (JMAX .LE. 191) GO TO 8 


aaaaaaaaan 
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JMAX = 181 
VTOP=VBOT + 188.*DELTV 
8 VVTOS=VTOP 


vo= VBOT 
VLWST=VBOT-VSTIP 
IP(VO.GT.2369.4) VLWST=2259.9 
VEGHST=VTOP + VSTEP 
DO 12 J=1,JIMAX 
1g OMEGB(I)=VO+ (J-1)*DELTV 
Cxe** (SUML,SUA2)=S(LEFT,RIGHT) SIDE OF REGION 
C ***Ze5RO SUM?S -- READ TAPE -- CALCULATE 
DO 12 MM=!,MAX 
DQ 12 NPT=1,NPIPTS 
DO 12 J=1,J0MAxX 
SUML(MM, NET ,J)=9.0 
SUM2 ("M,NPT,J)=9.9 


nN 


ar 


3 DO 3.5 MM=1,MAX 
# M=MS>EC (Mit) 
IF (M.LE.@' GO TO 38 
c *** LOGIC IS FOR TEMP DEP OF LINESTRENGTH 
LOGIC=.FALSE. 
TF(M.FQ.1.0R.M.EQ.3.0R.M.EQ.6) LOGIC=.TRUE. 
C*** PRESSJRE, TEMPERATURE LOOP 
pO 23 NPT=1,NPTPTS 
c*** (D,T) COMPUTE DEPENDENCE OF LINE PARAMETERS 
CSL (NPT) =(TO-T(NPT) )/(TO*T(NPT) *9.6946) 
WT=SQRT(TU/TCNPT) ) 
CS2 C1,NPT)=T@/T(MPT) 
TF(LCGIC) CS2(M,NPT)=CS2(M,NPT) *WT 
CA(N?T) =WI"P (NPT) /PY 
22 CONTINUE 
38 CONTINUE 


CALL REDLIN(VBOT,VTOP) 


Cc 
ses CALCULATION OF FAR LINES IS COMPLETED 

it=1 

DO 33 J=1,JMAX 

WRITE(6,34) SPECIE(M) ,OMEGB(J),(SUM1(M,NPT,J) ,SUM2(M,NPT,J) ,NPT#1 

a) 
34 FORMAT(AS,F19.3,5(2X%,2E9.3)) 

MeM+1 

TR(M.GT.5) M=5 
33 CONTINUE 
Cc 


C*** CALCULATE NEAR LINES 
c *** CALL STRONG FOR CALCULATIONAL POINTS (OMEGA?S) 
REWIND 13 
VTOP =V1+VBLOCK 
DO 355 IFILE=1,NFIL&S 
RBWIMID 11 
VBOT=V1-BOUND 
LE(VIOP.GT.VVTOP) VTOP=VVTOP 
IF(VTOP.GT.V2) VTOP=V2 
NREC=3 


tr 


348 CONTINUE 
C*cpC READ(13) N,(GNU(I),S(I),ALPHA(I),EDP(I),MOL(I), 
READ (13,END=359) N, (GNU(I),S(I),ALPHA(I),EDP(I) ,MOL(I), 
X LCHK(I),I=1,N) 
C*cpc IF(EOF(13)) 358,345 
345  WREC=NREC+1 
WRITE(1LL)N, (GNU(I),S(I) ,ALPHA(I),EDP(I),MOL(I), 
X LCHK(I),I=#1,N) 
GO TO 349 
35@ WRITE(4,382) V1,VTOP 
382 FORMAT(2F1@.2) 
CALL STRONG(MSPEC,V1,VTOP,DV,DV2,JCALC,SSTR,MAX,NREC) 
VTOP = VTOP+BOUND 
C*** ZERO STOR---READ TAPE --- CALCULATE 
DO 35 MM=1,MAX 
DO 35 NPT=1,NPTPTS 
JMAX=JCALC (MM) 
DO 35 J=1,JMAX 
35 STOR(MM,NPT,J)=9.2 
REWIND 11 
DO 55 IREC=1,NREC 
READ(11)N,(GNU(I),S(I),ALPHA(I),EDP(I),MOL(I), 
X LCHK(I),I=1,N) ; 
bU 48 NPT=1,NPTPTS 
DO 46 L=1,N 
M=MOL(L) 
DO 37 MM=1,MAX 
IF(M.EQ.MSPEC(MM)) GO TO 52 
37 CONTINUE 
GO TO 46 
52 SO=S(L)*CS2(M,NPT) *EXP(-EDP(L) *CS1 (NPT) ) 
AL=ALPHA(L) *CA(NPT) 
AD=.3581E-6*GNU(L)*WT 
ARAT= .83255*AL/AD 
c*** LOOP OVER CALCULATIONAL POINTS (OMEGA?S) 
K=1 
JMAX=JCALC (MM) 
DO 44 J=1,JMAX 
V=OMEGA(J,MM) 
Cc #** DETERMINE RANGE TOR ACCEPTING LINES 
38 VLEFT=OMEGB (K) 
IF (VLEFT+DELTV-V) 39,489,498 
39 K=K+l 
GO TO 38 
42 21=A35 (VLEFT-GNU(L) ) 
IF (21.LT.BOUND) GO TO 42 
Z2*ABS (VLEFT+DELTV-GNU(L) ) 
IF (22.GT.BOUND) GO TO 44 
c*** IF LINE WITHIN "BOUND" OF INTERVAL, INCLUDE IT IN CALCULATION 
42 Z=ABS (V-GNU(L) ) 
ETA=.83255*2Z/AD 
CALL ABSORB(AL,AD,ETA,ARAT,2Z,AK) 
FAC=SO*AK 
c** CO2 FORM FACTOR 
IF(M.EQ.2) FAC=FAC*FORM(2Z,M) 
STOR(MM,NPT,J)=STOR(MM,NPT,J)+FAC 


oe a 


. 
“ 


63 
64 


66 
68 
c 


CONTINUE 
CONTINUE 
CONTINUE 

MORE LINES 2??? 
CONTINUE 


. NEARBY REGIONS ARE NOW CALCULATED 


COMBINE RESULTS 
DO 89 MM=1,MAX 
M = MSPDEC(MM) 
K=t 
JMAX = JCALC (MM) 
DO 63 J=1,JMAX 
V=OMEGA (J , MM) 
VLEFT=OMEGB (K) 
IF (VLEFT+DELTV-V) 63,64,64 
K=K+l 
GO TO 62 
CONTINUE 
DO 66 NPT=1,NPiTPTS 
FAC1=SUM1 (MM, NPT,K) 
PAC2=SUM2 (MM,NPT,K) 


STOR(MM,NPT,J)=STOR(MM,NPT,J) +FAC1+(FAC2-FAC1)*(V-VLEFT) /DELTV 


STOR‘'MM,NPT,J) = STOR(MM,NPT,J) *W(M) 
CONTINUE 


CRARS ERY WRITE TABLE FOR MRDA 


c 


Cc 
C2 


22 


228 


aaa 


365 


WRITE (4,228) SPECIE(M) ,JMAX 
WRITE(6,222) SPECIE(M) ,JMAX 
FORMAT (1X,A4,15) 
FORMAT (A4,15) 
DO 74 J=1,JMAX 
WRITE(6,226) OMEGA(J,MM), (STOR(MM,NPT,J) ,.NPT=1,NPTPTS) 
WRITS (4,226) OMEGA(J,MM), (STOR(MM,NPT,J),NPT=1,NPTPTS) 
PORMAt (F12.3,9E12.6) 
CONTINUE 
CONTINUE 
V1i=V1+VBLOCK 
IF(IFILE.EQ.NFILES) Vi=VVTOP 
IF(V1.GE.V2) STOP 23 
VTOP=V1 + VBLOCK 
CONTINUE 
GO TO 5 


END 
FUNCTION FORM(2,M) 


FORM FACTOR FOR SUB-LORENTZIAN TAILS 


FOPM = 1.8 

If (M.NE.2) RETURN 

IF (2.LT.9.5) RETURN 
IF (42.6T.23.) SO TO 19 
FORM? 1.869 *EXP(-.133%*2) 
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4. 


18 


18 


15 
28 


c*CLC 
82 


84 
c*CDC 


81 
83 


85 


87 


809 
89 


RETURN 


FORM = .@5 ‘ 
IF (2.LE.5@.) RETURN 
FORM = 9.9 


IF (2.GE.258.) RETURN 

PORM= .@9@5*(12.5-.85*2) 

RETURN 

END 

SUBROUTINE ABSCRB(AL,AD,ETA,ARAT,Z,AK) 

IF ((ETA.LE.5.).AND.(ARAT.LE.2.)}) GO TO 19 
AK=(.31831)*AL/(Z**2 + AL*AL) 

GO TO 28 

CONTINUE 

AK=8.0 

DO 15 K=e1,51 

Maq20 4 (kate ot 
FY=(ARAT/SD)*.14952*EXP(-Y*Y)/(ARAT**2 + (ETA-Y)**2) 
PY=FY*.1 

AK=AK + FY 

CONTT.NUE 

RETURN 

END 

SUBROUTINE DATIN 

COMMON/INPUT/ P(18),T(18) ,W(7),V1,V2,DV,VLWST,VHGHST,DELTV 
K,BOUND,NPTPTS ,MSPEC(7),SSTR,VBLOCK,DV2 

READ(5,76) NPTPTS,MSPEC 

READ(5,76,END=345) NPTPTS ,MSPEC 

FORMAT (812) 

READ(5,77) (P(I),I=1,NPTPTS) 

READ(5,77,END=345) (P(I),I=1,NPTPTS) 

FORMAT (8(E18.@)) 

READ(5,77,END=345) (T(I),I=1,NPTPTS) 

READ(5,77) (T(I),I=1,NPTPTS) 

WRITE(6,82) (P(I),I=1,NPTPTS) 

PORMAT(* PRESSURE=*,5(2X,F7.2)/108X,5(2X,F7.2)) 
WRITH(6,84) (T(I),I=1,NPTPTS) 

FORMAT(* TEMPERATURE=*,5(2X,F7.2)/13X,5(2X,F7.2)) 
IF(EOF(5).NE.@) GO TO 345 

READ(5,81) (W(M),M=1,7) 

FORMAT (7E19.3) 

WRITE(6,83) 

FORMAT (3X, *WATER* ,6X,*CO2*,6X, *OZONE*, 7X, *N20*,7X,*CO*,8X, *CH4*, 
17X,*02*,4K%) 

WRITE(6,81) (W(M) ,M=1,7) 

R5AD(5,85) V1,V2,DV,VLWST,VHGHST,DELTV,BOUND 

FORMAT (7518.3) 

WRITE‘6,87) Vi,V2,DY,VLWST,VHGHST,DELTV,BOUND 
FORMAT (* (V1,V2,DV) =*,3F18.3,5X,* (VLWST,VHGHST,DELTV) =#=*, 
K3F1@.3,5X,*BOUND =*,F18.3) 

READ(5,889) SSTR,VBLOCK,DV2 

FORMAT (E1J@.2,F18.3,3F18.2) 

WRITE(6,89) SSTR,VBLOCK ,DV2 

FORMAT (* SMIN FOR STRONG = *,E1@.3,5X,*VBLOCK =*,F1@.1,5X, 
&*DV2 =*,F5.2) 

TF (V1.GE.V2) STOP 21 ‘ 
RETURN 
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345 STOP 28 
END 
SUBROUTINE STRONG (MSPEC,V1,V2,DV,DV2,JCALC,SSTR,MAX ,NREC) 
C DETERMINES CLACULATIONAL POINTS FOR STRONG LINES 
LOGICAL LCHK,LCHKS 
DIMENSION G@(281),S@(281).,LCHK#(2@1) ,DENS(7) ,MSPEC(7) 
DIMENSION JCALC(7) ,SPECIE(7) ,M@(281) ,A®(281) 
f COMMON/LINES/ GNU(25@),S(258) ,ALPHA(258) , EDP(258) ,MOL(25#), 
; KLCHK(25@) ,TI(258,9),ITI(250) , TMAX 
; COMMON/OMG/ OMEGA(291,6),STOR(6,9,384) 
DATA DENS/1.@,1.8,.81,.081,.8085,.085,.881/ 
DATA SPECIE/4HH20 ,4HCO2 ,4HO3 ,4HN20 ,4HCO ,4HCH4 ,4HO2 / 
c*** INITIALIZE CONTROL VARIABLES 
W18 = SORT(19.9) 


NMAX = 288 
c** TOO MANY LTNES---~INCREASE SMIN AND READ TAPE AGAIN 
5: NLINES=9 

REWIND 11 


DO 3@ IREC=1,NREC 
; READ (11)N, (GNU(I),S(I) ,ALPHA(I) ,EDP(I) ,MOL(I), 
X LCHK(I) ,I=1,N) 
IF(GNU(N).LT.V1) GO TO 39 
C** THROW OUT LINES THAT ARE WEAKER THAN SSTR 
DO 20 J=1,N 
IF(GNU(J).LT.V1) GO TO 29 
IF(GNU(J).GT.V2) GO TO 27 
IF (S(J).LT.SSTR .AND. (.NOT.LCHK(J))) GO TO 28 
NLINES = NLINES+1 
IF(NLINES.LE.NMAX) GO TO 15 
SSTR=W10*SSTR 
GO TO 5 
15 G@(NLINES) =GNU(J) 
$0 (NLINES) =S (J) 
Ag (NLINES ) =ALPHA (J) 
M@ (NLINES ) =MOL‘J) 
LCHK9 (NLINES ) =LCHK (J) 
2@ = CONTINUE | 
38 CONTINUE 
Cc *** SET UP PARAMETERS 
27 ID@=5 
DELT1 = ID@*DV 
ISPACE = 2*IDZ + 2 
LMAX=2088 
LMAX1=LMAX - ISPACE - 5 
LMIN= LMAX/2 
¢ *** PICK OUT STRONG LINES 
FAC=1.9E-28 
JLOOP=9 
389 JLOOP=JLOOP+1 
IPTS=0 
FAC=.3*FAC 
DO 37M M=1,MAX 
JMAX = LMAX-JPTS 
JCALC (M) = 
IF (MSPEC(M).LE.@) GO TO 378 
OMEGA(1,M)#V1 
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MM=MS PEC (M) 
SMIN=FAC*DENS (MM) 
JCNT=1 
DO 328 JC=1,NLINES 
IF(M@(JC).NE.MM) GO TO 328 
c*** GNU.GT. 2388 --- ACCEPT ALL CO2 R BRANCH FUNDAMENTAL LINES 
IF (G#(JC).GT.2388..AND.LCHK@(JC))GO TO 388 
IF(S9(JC).LT.SMIN) GO TO 32g 
c*** ACCEPT LINE - DETERMINE CALCULATIONAL POINTS 
388 D1l=(G@(JIC)-OMEGA(JCNT,M))/DV 
IF (D1.LT.(.5*DV)) GO TO 328 
ID=IDg 
IF (D1.GE.ISPACE) GO TO 313 
ID=(D1+@.5) 
c*** ID IS NUMBER OF POINTS ON LINE WING---DV=SPACING BTWN PTS 
IF (ID.EQ.@) ID=1! 
D2=(CH(IC)-OMEGA(JCNT,M))/ID 
JCNT=JCNT+ID 
OMEGA (JCNT,M)=G8 (JIC) 
c*** IF (JCNT.GT.190) JLOOP=69 
IF ({D.EQ.1) GO TO 328 
ID=ID-1 
GO TO 316 “ 
c*** ID.GE.8---LINES WELL SEPARATED - L.H.S OF INTERVAL 
313 DO 314 II=1,ID 
JCNTII=JCNT+II 
314 OMEGA(JCNTII ,M) =OMEGA(JCNT,M)+II*DV 
C*** DETERMINE POINTS BETWEEN THE TWO LINES---INCREMENTS OF Dv2 
JCNTID=JCNT+ID 
VA = DV2*(FLOAT(INT(OMEGA(JCNTID,M)/DV2)) + 1.8 ) 
VB = DV2*FLOAT(INT( (G@(JC)-DELT1-.@81)/DV2)) 
JCNT = JCNT+2*ID+1 
IF (VB.LT.VA) GO TO 316 
ID2 = 1+INT((VB-VA)/DV2+.9995) 
JCNT = JCNT-ID-1 
DO 315 II=1,ID2 
JCNTII=JCNT+II 
315 OMEGA(JCNTII,M) = VA+(II-1)*Dv2 
JCNT = JCNT+ID+ID2+1 
316 OMEGA (JICNT ,M) =GH(IC) 
c*** LINES CLOSE --- OR ~-- R.H.S. OF INTERVAL 
DO 317 II=1,ID 
JCNTIL=JCNT-II 
317 OMEGA(JCNTII,M)=OMEGA(JCNT,M)-II*DV 
Cane IF TOO MANY CALC.PTS., INCREASE SMIN 
IF (JCNT .GT. JMAX) GO TO 371 
328 CONTINUE 
c*** SCALC = NO,. OF CALCULATIONAL POINTS 
JCALC (M) =JCNT 
378 JPTS=JPTS+JICALC (M) 
Grr DECIDE WHETHER TO LOOP BACK AGAIN 
‘IF (J7,00P.GE.6@) GO TO 371 
IF ((JPTS.LT.LMIN) .AND.(FAC.GE.SSTR)) GO TO 389 
IF (JPTS.LT.LMAX1) GO TO 321 
371 FAC=4.*FAC 
GO TO 309 
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321 CONTINUE 
ICOUNT=% 
DO 248 M=1,MAX 
MM=M5 PEC (1) 
IF (MM.LE.@) GO TO 248 
JCNT=JCALC (M) 
IF(JCNT.GT.1) GO TO 338 


NO STRONG LINES IN THIS BLOCK 


JCALC(M)=5 
JCNT = 5 
DLT = .25*(V2-V1) 
DO 331 K=2,5 
331 OMEGA(K,M)=V1+(K-1)*DLT 
WRITE(6,339) JLOOP,FAC,SPECIE(MM) , (OMEGA(JJ,M),JJ*1,JCNT) 
339 FORMAT(/* NO INTENSE LINES*,106X,*JLOOP =*,I4, 
&10X,*SMIN=*,E19.3,12X,*SPECIES\* ,A1@/ 
&* OMEGA=*,5F14.3) 
GO TO 248 
Cc 
c STRONG LINES IN THIS BLOCK 
c*** DEFINE CALCULATIONAL POINTS NEAR V2 
338 Di=(V2-OMEGA(JCNT,M))/DV 
c** LAST LINE VERY CLOSE TO V2???? 
IF (D1.GT.1.8) GO TO 341 
JCNT = JCNT+1 
OMEGA(JCNT,M) = V2 
GO TO 335 
ID=1ID8 
IF (D1.GE.ISPACE) GO TO 343 
ID=(D1+H.5) 
IF (ID.EQ.@) ID=1 
D2=(V2-OMEGA(JCNT,M))/ID 
JCNT=JCNT+ID 
OMEGA (JCNT ,M)=V2 
IF (ID.EQ.1) GO TO 335 
ID=ID-1 
GO TO 346 
DO 344 II=1,ID 
JCNTIL=JCNT+II 
OMEGA (JCNTII,M) =OMBGA(JCNT,M)+II*DV 
SCNTINSICNT+ID+1 
OMEGA (JCNTID,M) =(OMEGA(JCNT,M)+V2)/2. 
JCNT=JCNT+2* (ID+1) 
OMEGA (JCNT,M)=V2 
346 DO 347 II=1,ID 
JCNTII=JCNT-II 
347 OMEGA(JCNTII,M)=OMEGA(JCNT,M)-II*DV 
C*** PRINT STRONG LIVES AND CALC.PTS. 
335 JCALC (M) =JCNT 
WRITE(6,325) JLOOP,JCNT,FAC,SPECIE(MM) 
325 FORMAT (/* MOST INTESE LINES:*,1@X,*JLOOP =*,I3,5X, 


*I5,* CALCULATIONAL POINTS*,5X, 
**SMIN=*,E18.3,18X,A18/5X,*V*,10X,*S*,6X, *ALPHA* , 3X, *M* ,18X, 
**INTERMEDIATE POINTS*) 

JCMIN=1 

JMAX1=JCNT-1 

J2 = ISPACE-1 

WRITE(6,327) (OMEGA(J,M) ,J=1,ISPACE) 

FORMAT (F10.3,26X,11F9.3,(/35X,18F9.3)) 

DO 368 JJ=2,JMAX1 

DO 328 JC=JCMIN,NLINES 

IF (OMEGA(JJ,M).NE.G@(JC)) GO TO 328 

Jil=JT+I2-1 

WRITE(6,326) GO(JC),SH8(IC) ,AB(IC) ,MO(JIC), (OMEGA(S+1,M),J=J7,71) 
JCMIN = JC 

GO TO 329 

CONTINUE 

CONTINUE 

CONTINUE 

FORMAT (F1%.3,2E9.2,13,5X,11F9.3,(/35X,18F9.3)) 
CONTINUE 

RETURN 

END 

SUBROUTINE REDLIN(Vi,V2) 

DIMENSION CS1(9),CS2(7,9),CA(9) ,SUM1(6,9,1M1) ,SUM2(6,9,1@1) 
DIMENSION OMEGB(191) 

DIMENSION TI(25@,9) ,ITI(25@) 

DIMENSION GNU(25@) ,S(25@) ,ALPHA(25@) ,MOL(25@) , LCHK (258) 
DIMENSION EDP(25@) ,DENS(7) ,CO2R(5) ,MSPEC(7) ,ALPHB(6) 
COMMON/LINES/GNU,S,ALPHA,EDP,MOL,LCHK,TI,ITI,TMAX 
COMMON/INPUT/P(18) ,T(18) ,W(7) ,VA,VB,DV,VLWST,VHGHST,DELTV, 
X BOUND,NPTPTS ,MSPEC,SSTR,VBLOCK ,DV2 
COMMCN/BLOCK1/MAX,CS1,CS2,CA,SUM1 ,SUM2 ,OMEGB ,JMAX ,NFILES 
LOGICAL TOEND,IOND13 ,LCHK 

DATA DENS/1.8,1.0,8.01,8.8981,8.8805,8.885,8.9081/ 

DATA CO2R/6H ® 8,6H @ 1 1,6H 6538 8B ,6HS 1 
DATA ALPHB/9.87,8.11,8.08,8.86,8.855,8.848/ 

IFILE=9 

IOEND=.FALSE. 

IOND13=.FALSE. 

REWIND 11 

END FILE 11 

IN=11 

IOUT=12 

REWIND 2 

REWIND IN 

REWIND IOUT 

REWIND 13 

VLOW=V1-BCUND 

VHISH=V2+BOUND 

VL=Vi-VBLOCK 

VL=ViL+VBLOCK 

IF(VL.GE.V2) GO TO 238 

VR=VLL+VBLOCK 

IF(VR.GT.V2) VR=V2 

VBOT=VL-BOUND 

VTOP=VR+BOUND 


125 


999 


138 
C*cDC 


c*cDc 
158 


VRB=VR-BOUND 

tel 

IF(I.LE.25@) GO TO 138 

N=250 

WRITH(6,999) VL,VR,VBOT,VTOP,VRB,VLWST,VHGHST,VSTEP,VLOW, VHIGH 
FORMAT (1X,1@F19.3) 

WRITE(13)N,(GNU(I),S5(I),ALPHA(I),EDP(I) ,MOL(I),LCHK(I),I=1,N) 
I=1 

CONTINUE 

READ(IN)GNU(I),S(I),ALPHA(I),EDP(I) ,MOL(TI) ,LCHK(T) 
READ(IN,END=219)GNU(I),S(I),ALPHA(I),EDP(I) ,MOL(I) ,LCHK(TI) 
IF(EOF(IN)) 217,158 

V=GNY(T) 

IF(VR.GE.V2) GO TO 169 
IF(V.GE.VRB)WRITE(IOUT)GNU(I),S(I),ALPHA(I) ,EDP(I) ,MOL(TI), 


xX LCHK (I) 


168 


| 
| 
' 
| % LCHK(I),I=1,N) 


288 
218 
215 


228 


238 
c*cpCc 


c*cpc 


IF(V.LE.VTOP) GO TO 280 

TF(IOEND) GO TO 139 

N=I-1 

WRITE(6,999) VL,VR,VBOT,VTOP,VRB,VLWST, VHGHST,VSTEP , VLOW, VHIGH 
IF(N.GT.@)WRITE(13)N, (GNU(I),S(I) ,ALPHA(I) ,EDP(I) ,MOL(I), 


END FILE 13 

IFILE=IPILE+1 

IOEND=.TRUE. 

GO TO 138 

T=I+1 

GO TO 125 

IF(IOEND) GO TO 228 

IF(.NOT.IOND13) GO TO 239 

IF(VR.LT.V2) GO TO 238 

NeI-1 

WRITE(6,999) VL,VR,VBOT,VTOP,VRB,VLWST,VHGHST,VSTEP,VLOW, VHIGH 
IF(N.GT.@)WRITE(13)N, (GNU(I),S(I),ALPHA(I),EDP(I),MOL(I), 
X LCHK(I),i=1,N) 

END FILE 13 

IFILE=IFILE+1 

GO TO 238 

END FILE I0UT 

ISAVE=IOUT 

ITOUT=IN 

IN=ISAVE 

REWIND IN 

REWIND IOUT 

IOEBND=.FALSE. 

GO TO ine 

CONTINUE 

READ(2) TMIN,TMAX,NREC, ((TI(I1,031),031°1,9) ,DUM, 
READ(2,END=28@,ERR=248) TMIN,TMAX,NREC, ((TI(1I1,d1),01*1,9),DUM, 
X DUM,DUM,ITI(I1),1I1=1,NREC) 

JTO=IOCHEC (2) 


C*CDC IF(JTO) 249,278 


278 


CONTINUE 


C*CDC IF(EOF(2)) 288,298 
GO TO 299 
248 WRITE(6,245) GNU(I) 


245 FORMAT(* PARITY ERROR ENCOUNTERED AT*,F1i2.3) 


GO TO 238 


288 IEOF=IROF+1 


NEOF=NEOF+1 
IF(NEOF.GT.2) STOP 22 
GO TO 238 


298 N=EOF=8 


IF(TMAX.LT.VLWST) GO TO 238 
IF(TMIN.GT.VHGHST) GO TO 1588 

DO 119@ K=1,NREC 

DO 319 MM=1,MAX 
IF(ITI(K).EQ.MSPEC(MM)) GO TO 328 


318 CONTINUE 


GO TO 1198 


328 V=TI(K,1) 


IF(V.GT.VHGHST) GO TO 1118 

IF(V.LT.VLWST) GO TO 1188 

M=ITI(K) 

SMIN=TI(K,2)*DENS (M) 

IF(SMIN.LE.1.9E-27) GO TO 1198 

IF(I.LE.25@) GO TO 388 

IF(IOEND) GO TO 388 

N=259 

WPRITE(6,999) VL,VR,VBOT,VTOP,VRB,VLWST,VHGHST,VSTEP,VLOW, VHIGH 
WRITE(13)N, (GNU(I),S(I),ALPHA(I) ,EDP(I) ,MOL(I),LCHK(I), 
X I=1,N) 

T=1 


388 GNUCI)=V 


385 


S(I)=TI(K,2) 

ALPHA(I)=TI(K,3) 

IF(M.£Q.1) GO TO 398 

IF(ALPHA(I).GT.@.8) GO TO 385 

ALPHA (I) =ALPHB(M-1) 

GO TO 398 

IF(ALPHA(I).LT.#.@1.OR.ALPHA(I).GT.1.@) ALPHA(I)=8.86 


398 EDP(I)=TI(K,4) 


MOL(I)=M 

IF(M.NE.2) GO TO 1888 
LCHK(I)=.FALSE. 

DO 499 J=1,5 
IF(CO2R(J).NE.TI(K,J3+4)) GO TO 1888 


488 CONTINUE 


LCHK(I)=.TRUE. 


1888 IF(V.LT.VLOW) GO TO 1855 


IF(IOND13) GO TO 1@55 

IF(VR.GE.V2) GO TO 1828 

IF(V.GT.VHIGH) GO TO 182@ 

IF(V.GE.VRB) WRITE(IOUT)GNU(I),S(I),ALPHA(I),EDP(I) ,.MOL(I), 
X LCHK(T) 


1828 IF(V.LE.VTOP) GO TO 1858 


IF(V.GT.VHIGH) IOND13=.TRUE. 

IFC(IOEND) GO TO 1855 

WRITE(6,999) VL,VR,VBOT,VTOP,VRB,VLWST,VHGHST,VSTEP,VLOW, VHIGH 
NsI-t 

IF(N.GT.@) WRITE(13)N, (GNU(I),S(I),ALPHA(I),EDP(I),MOL(I), 
X LCHK(I),I#1,N) 
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1858 


1855 
1857 


1588 


END FILE 13 


IFILE=IFILE+1 
IOEND=.TRUE. 
GO TO 1955 
L=I 

I=I+1 

GO TO 1957 
LeI 


DO 198@ NPT=1,NPTPTS 
$O0=S(L) *CS2(M,NPT) *EXP(-EDP(L)*CS1 (NPT) ) 
AL=ALPHA(L) *CA(NPT) 

DO 1978 J=i,JMAX 
VV=OMEGB (J) 

Z1=ABS(VV-V) 
IF(Z1.LT.BOUND) GO TO 187% 
22=ABS (VV+DELTV-V) 
IF(Z2.LT.BOUND) GO TO 197% 
FRM1=1.9 

FRM2=1.8 

IF(M.NE.2) GO TO 1869 
FRM1=FORM(21,M) 
FRM2=FORM(22 ,M) 

CONTINUE 

SUM1 (MM,NPT,J)=SUM1(MM,NPT,J)+FRM1*.3183*SO*AL/ (Z1**2+AL**2) 
SUM2 (MM,NPT,J)=SUM2 (MM,NPT,J) +FRM2*.3183*SO*AL/ (Z2**2+AL**2) 
CONTINUE 

CONTINUE 

CONTINUE 

IF(.NOT.IOEND) GO TO 239 
END FILE IOUT 

ISAVE=IOUT 

IOUT=IN 

IN=ISAVE 

IOEND=.FALSE. 

REWIND IN 

REWIND IOUT 

GO TO 198 

NFILES=IFILE 

RETURN 

END 
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APPENDIX D 
MIDTRAN FLOWCHART; SUBROUTINE LIST 


Table D-1, Generalized MIDTRAN Flowchart 


20 Read Data Files If Model < 0 Yes 
Model! Atmos phe res No 
180 Broadband Data 


Call LIB 


190 STOP 22 
Calculate Horizontal | 200 Trajectory 
at Profile for Each Layer Misses Earth 


P 
Calculate Vertical Profiles rue 


and Equivalent Sea Level 
680 Absorber Amounts 


690 


Broadband Calculations: 
Aerosol Extinction 
880 N2, H20 Continua 


920 
Print Layer Profiles For 
High Resolution Section 


sas High Resolution 
Calculations: Read Tape; 
Call PTPTS; 
Monochromatic Calculation 


i Combine High Resolution 
and Broadband Calculation 


1200 Store Results On 
Output File 9 


Table D-2, Listing of MIDTRAN Subroutines 


; 


Computes the equivalent absorber amount at a given 
altitude (from LOWTRAN). 


Determines the interpolation points from the data tape 
for each layer. 


Computes the zenith angle given H1, H2 and the earth 
center angle (from LOWTRAN). 


Performs output functions 
Reads Cards 5, 6, 7. 


Removes trailing blanks in Card 5. 


Degrades results to the desired resolution using the 
generalized slit function. 


Sets up frame for the plots 
Prints output and plots curves. 


Skips over data sets on File 9. 
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APPENDIX E 


LIST OF MIDTRAN VARIABLES 


CON 


Width of fixed triangular slit function ( = 0.1 em™2) 


Aerosol number density for MODEL = 1 
Aerosol number density for MODEL = 2 


Extinctica coefficient read from tape for Kth pressure- 
temperature point at frequency VV 


Interpolated extinction coefficient 

Equivalent absorber amount per km at level L 
Wavelength (um) 

Altitude at level Z (LYR) 

Input zenith angle (degrees) (compare with 6, in the text) 
Average wavelength used in refractive index expression 
Constant A defined in Eq. (10) of LOWTRAN3 Manual 
Blackbody function 


Angle subtended at the earth's center as path transverses 
adjacent levels 


Total angle subtended by path at earth's center 
Conversion factor from degrees to radians 


Wavelength dependent coefficient used in refractive 
index expression 


Species concentrations 


Wavelength dependent coefficient used in refractive 
index expression 


Absorption coefficient for nitrogen (~ 4 um) 
Absorption coefficient for water vapor continuum 
Extinction coefficient for molecular scattering 
Extinction coefficient for aerosol models 
Aerosol absorption coefficient 

Water vapor amount (pr.cm/km) at level L 
Increment for fixed slit function 

Optical depth of a species 

Dew point temperature (°c) 

Path length from level Lto Level L+ 1 


Wavenumber increment 
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IC 

ICNT 
ICOUNT 
IDV 
IFIND 


MRDA frequency interval 


Height increment from level L to level L + 1 
Equivalent absorber amounts per km at height H 


Equivalent absorber amounts 


Emissivity of background radiation source 
End-of-file control variable 
Integrated absorber amount from level Lto levelL + 1 


Dummy variable 


Summing variable for transmittance 
Log transmittance for radiation 
Summing variable for transmittance 


Frequency for print-out 


Intermediate result in interpolation of AK 
Degraded radiation 
Dummy frequency variable 
Intermediate result in interpolation of AK 
Altitude 

3 


Aerosol number density (no. cm ~) 


Estimated tangent height (km) 
Minimum altitude of path trajectory (km) 


SI ae ee eee ee NOY ena ae 


Dummy variable for transmittance 


Aerosol number density (no. em”) for 23 km 
visual range 


Aerosol number density (no. em”) for 5 km 
visual range 


Initial altitude (km) 

Final altitude (km) 

Counting variable for low resolution calculations 
Index for low resolution calculations 

Index for low resolution printed output 


Frequency increment 


Call parameter for subroutine ANGL 


Fg TE TT EE 


Aerosol model indicator 


Integer variable for printing heading 


Parameter used when reading in a new atmospheric a 
model 

Indicator for using subroutine POINT to calculate 

refractive index only (IP = 0) or equivalent absorber * 


amounts also (IP # 0) 

Radiation calculation flag 

Indicator for type of atmospheric path 

Frequency of calculations 

Starting frequency 

Last frequency 

Running integer for altitude identification 

Print option parameter 

Level indicator for altitude H1 

Level indicator for altitude H2 

General loop variable 

Elements in P-T matrix used for AK interpolatim 
Number of species for high resolution calculation 6) 


Cycling parameter for downward path 

Running index for layers 

Total number of levels transversed in the path 
Parameter used for defining longer of two paths 
Running index for level s 

Counting variable for long path storage 


Number of layers for low resolution radiance 
calculations 


Counting variable for layer index 
Altitude of Lth layer in path 
Frequency identifier for transmittance calculation 


Frequency identifier for transmittance calculation 


Index for model atmosphere 


Number of radiation and transmission calculations 
Integer used to identify required model atmosphere 


Variable to select temperature profile and 
counting variable 


Variable to select water vapor profile 
Variable to select ozone profile 

Loop variable 

Plot control parameter 

Plot control parameter 


Frequency indicator for water vapor continuum 
transmittance calculation 


Number of levels in model atmosphere data 
Number of layers in model atmosphere data 
Number of points in the pressure temperature matrix 
Determines units for radiation output 

Control variable for radiation slit function 
Determines scale for transmittance 

Control variable for transmittance slit function 
Index for locating species on library tape 
Pressure 

Angle of arrival at H2 

Angle of arrival at H2 

3. 141592654 that is (1°) 

Pressure values on library tapes 

Partial pressure H90 

Pressure at level LL 

Total pressure in atmospheres 

Angular deviation of path from initial direction 
Radiance 

Path length (km) 

Earth radius (km) 

Refractive index of air at level L 
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RH 
RN 


Relative humidity (%) 


Ratio of refractive indices of air above and 
below a given level 


Ratio of earth center distances between adjacent 
levels 


Sine of angle of arrival at adjacent level (cf sin a) 
Number of species 


Sine of the local zenith angle at a given level 
(cf sin 6) 


Slant range (km) 

Sum of optical thicknesses of absorbers 4 thru 8 
Accumulated integrated absorption 
Temperature (°K) at level L 

Transmittance 

Transmittance 

Background radiation calculation temperature 
Temperature at level LL 

Zenith angle at a given level (in radians) 
Zenith angle at a given level (in degrees) 
Ambient temperature ("Cc ) 

Total transmission 

Transmittance from fixed slit function 
Broadband transmittance of layer LL 

Ratio 273. 15/(TMP + 273. 15) 

Temperature values on library tape 

Total transmittance at frequency V 
Absorption due to aerosol only at frequency V 
Refractive index of layer above initial altitude H1 
Frequency (em™2) 

Initial frequency in tape data block 

Final frequency in tape data block 


Used to compare lower frequency of tape data block 
with calculation frequency 
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VCHK2 Used to compare upper frequency of tape data block 
with calculated frequency 


- VH Integral of the equivalent absorber amounts from 
H1 to level L 
VS Visual range (km) at sea level 
i . VMAX Max frequency contained in tape 
f VMIN Minimum frequency contained in tape | 
VT Frequency for fixed slit function 
VV Frequency array read from tape 
vvl1 Used in interpolating tape input frequencies to 
calculation frequency 
VV2 Used in interpolating tape input frequencies to 
calculation frequency | 
Vx Wavelength at which aerosol coefficients are | 
read in (um) 
Vo Initial calculation frequency 
Vil Initial frequency for transmittance calculation, em! 
V2 Final frequency for transmittance calculation, em! 
WwW Total equivalent absorber amount for entire path 
WGAS Gas concentration 
WH Water vapor density at level L (gm m7) 
WH20 Water vapor concentration | 
wo Ozone density at level L(gm m7) 
WO3 Ozone concentration 
WW Equivalent absorber amount from observer to 
level L 
Ww2 Water vapor density for atmospheric model M 
at level L + 1 (gm m7) 
Xx Input height to POINT subroutine 
XD Wavenumber interpolation parameter 
XH Wavenumber interpolation parameter in HO 


continuum calculation 
XI Wavenumber interpolation parameter | 
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X-coordinate for lower left corner of plot 
Aerosol extinction coefficient 

Earth center distance of level L 

Earth center distance of level L + 1 

Input zenith angle in radians 


Refractive index of layer below input height 
from POINT subroutine 


Y-coordinate for lower left corner of plot 


Aerosol absorption coefficient of frequency V 
Altitude at level L in km 
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